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DEFINITION OF SYMBOLS 


Symbol 

A 

A D 


M 


w 


m 


V. 

1 


V, 


w 


pw 


Definition 

propellant tank wall surface area 
pressurant gas distributor area 
dimensional constant 
pressurant molecular weight 
ambient temperature 
propellant temperature 
ullage mean temperature at cutoff 
pressurant inlet temperature 
propellant tank volume 
initial ullage volume 
propellant volumetric drain rate 
wall thickness 
pressurant specific heat 
wall specific heat 


g dimensional constant 

c 

h ambient heat transfer coefficient 

a 

k pressurant thermal conductivity 

p ullage pressure 

r propellant tank characteristic radius 

(maximum radius for cylindrical tanks) 


ft 2 

ft 2 

lb f ft/Btu 

lb /lb mol 
m 

°R 

°R 

°R 

°R 

ft 3 

ft 3 

ft 3 / sec 
ft 

Btu/lb °R 
m 

Btu/lb °R 
m 

lb ft 
m 

lb f sec 2 
Btu/sec ft 2 °R 

Btuft/sec ft 2 °R 

lb^/ft 2 

ft 


U= 


v 



DEFINITION OF SYMBOLS (Concluded) 


Symbol Definition 

R universal gas constant 

a gas compressibility 

p pressurant viscosity, M/L 

p wall density, M/L 3 

w 

6 time of pressurization 


lb f ft/°R lb mol 
(-) 

lb /ft sec 
m 

lb /ft 3 
m 

sec 


vi 


PREDICTION OF PROPELLANT TANK PRESSURIZATION 
REQUIREMENTS BY DIMENSIONAL ANALYSIS 

SUMMARY 


Pressurant gas requirements for launch and space vehicles may be predicted by 
analytical models of the pressurization process. However, preliminary design studies 
require a fast and reasonably accurate method of predicting without resorting to com- 
puter programs. Therefore, dimensional analysis of a large number of pressurization 
tests and computer runs was applied to develop an equation that predicts pressurant 
requirements for cylindrical and spheroidal propellant tanks with an accuracy of ±10 per- 
cent. 


INTRODUCTION 


Although the most accurate method of predicting pressurant requirements is with a 
computer program that has been verified by experiments, it is advantageous to have a 
fast, reasonably accurate method to determine the total mass of pressurant gas required 
without resorting to the computer. This type of analysis is necessary in comparison and 
optimization studies for preliminary design where the number of possibilities to be con- 
sidered precludes a detailed computer analysis of each case. This report presents a 
single, general expression for the total required mass of pressurant. The expression 
was developed by dimensional analysis of the results of about 30 pressurization tests and 
120 simulations on an IBM 7094 Computer. 


DIMENSIONAL ANALYSIS OF TEST AND COMPUTER RESULTS 


The total mass of pressurant gas required is a function of the ullage mean tem- 
perature at cutoff derived with the gas equation of state, 

pVM 

W Total a RT (1) 

m 

Therefore, the total pressurant mass required may be calculated if the ullage mean tem- 
perature at cutoff can be determined. In the most general case, the ullage mean 
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temperature at cutoff will be a function of 12 system design variables, seven physical 
properties, the mechanical equivalent of heat, and the gravitational constant: 


T m ” f(J ’ g c’ M w’ k ’ < y r ’ T o- T L- V V ’ A ’ h a’ T a’ °pw’ p w' 


d w' ■>- V l> V L’ V 


The 19 variables can be expressed in six fundamental dimensions, length (L) , 
mass (M), time (0), temperature (T), heat (H), and force (F). The two dimensional 
constants, J and g c , are included because heat and force can be expressed in terms of 
the other four dimensions. The dimensions of each variable are given in Table IV. 

Since any equation representing physical phenomena must be dimensionally homo- 
geneous, it must be possible to write equation (2) in a nondimensional form. Therefore, 
using 7 r as a symbol for a dimensionless group, equation (2) may be written as follows: 

7 Ti = f (7T 2 , 7 T 3 7T ) (3) 

where the group 7 q contains the ullage mean temperature at cutoff. 

According to Buckingham's theorem, the number of dimfensionless groups in 
equation (3) is given by: 

i = n - m (4) 

where n = number of variables 

m = maximum number of these variables which will not form a dimensionless 

% 

group 

The quantity m is often equal to the number of fundamental dimensions. The total 
number of variables in equation (2) , including J and g c , is 22, and the number of funda- 
mental dimensions is six. By trial and error, it was determined that the maximum 
number of variables in this case that will not form a dimensionless group is indeed six, 
the number of fundamental dimensions. Therefore, according to equation (4) there will 
be 16 dimensionless groups in equation (3). 

Dimensionless Groups 

Because the maximum number of variables that will not form a dimensionless 
group is six, the 16 dimensionless groups may be determined by choosing six variables 
to be common to all groups and, in turn, adding each of the remaining 16 variables to the 
first six to form 16 dimensionless groups. Although any six variables could be chosen as 
the six to be common to all groups, intuition leads to the choice of the two dimensional 
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constants (J and g c ) , the pressurant molecular weight (M w ) , thermal conductivity (k). 
and viscosity (p) and the propellant tank characteristic radius (r) . Adding the ullage 
mean temperature at cutoff (T m ) to the six common variables yields the first group: 


T a b ,,c , d e f _g 
ii = J g M k p r T 
1 c w m 


( 5 ) 


Now, substituting the dimensions: 

»br “ic 


n rFLTrMLTr ithITmITxTLT 

Lrra L J b^J L^J L J L J 


( 6 ) 


Since the quantity on the left of equation (6) is dimensionless, the exponent of 
each of the six dimensions on the right must be zero. This condition determines six 
simultaneous equations, 


[L] 

0 = 

a + b - d 

[M] 

0 = 

b + c + e 

[0] 

0 = 

-2b - d - 

[T] 

0 = 

-d + g 

[H] 

0 = 

-a + d 

[F] 

0 = 

a - b 


(?) 


Simultaneous solution of these equations yields the following: 


ing: 


a = g 

d= g 

b = g 

e = -3 g 

c = 2 g 

II 

1 

Oq 

Because the value of g is 

a = 1 

d = 1 

b = 1 

e = -3 

c = 2 

f = -4 
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Substituting these values into equation (5) yields the first group 


J g M 2 k T 
c w m 


The second group 7r 2 is found in the same manner by adding the pressurant 
specific heat C p to the six common variables. Each of the remaining variables is added 
one at a time to the six common variables, and the remaining dimensionless groups are 
determined in the same manner as the first two. Then, 


c n 



(9) 

J g M 2 k T 

c W 0 

n - r i 

(10) 

J g M* k T, 

n ~ „» r i 

(11) 

f‘ re T 
n5 = ~M ' 

(12) 


w 


V 


lco 

1 

II 

(13) 

A 


~ r 2 

(14) 

! 

1 

o 

tJ 

(15) 

8 k 

P r 3 
w 

(16) 

71-9 “ M 

w 

d 

w 

^lO — r 
V. 

l 

*11 ~~p~ 

(17) 


(18) 


4 



(19) 


* 12 “ 


g M p 
c w 


M 2 r 


- 


h r 
a 


^14 - 


Jg M 2 k T 
c w a 

p 3 r 4 


M V T 
w L 

7715 “ n r 4 


D 


*16 - r 2 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


Now, the groups ir 1# 7t 4 , and 7r 14 may be divided by 7r 3 with no loss of generality, 
because none of the 16 groups is eliminated in the process. Thus 7 r 4 , 7r 4 , and 7r 44 may 
be replaced 7 r^, and 7rJ 4 , respectively ,where 



In the same manner, the group 7r 11 may be replaced with n' n , where 
g M P V. 

t c w 1 , nA , 

* n = ^il ^12 = ( 24) 

This change occurs because the initial ullage mean temperature is proportional to the 
product pVj. 

It is possible to reduce the number of dimensionless groups necessary in this 
case by realizing that the wall specific heat, density, and thickness can enter the prob- 
lem only in the combination Cp W p w d w which is the wall heat capacity per unit area. 
Therefore, the groups 7 t 8 , 7t 9 , and 7r 10 may be combined into a single group 7r'. 
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(25) 


*Z “ *8 *9 *10 ~ 


l~i r‘ 


M k 
w 


(C p d ) 
pw w w 


Similarly, it is possible to combine 7r lg and ir 16 . Thus, 




‘15 


7T. 


•16 


M V T 
w L 

V r! 


(26) 


Further, the propellant tank volume and wall surface area should enter the prob- 
lem as the area-to-volume ratio (A/V) . Therefore, the group n 6 and n 7 may be com- 
bined into a single group n. 

6 


t 





(27) 


But since the approximate area-to-volume ratio is given by 


A ^ 2jt r 1 _ 2 

V ~ 7r r 2 i r (28) 

the group tt[ is nearly a constant and may be eliminated. 

D 

Finally, since the Prandtl number (ir 2 ) does not vary greatly, this group may also 
be dropped. Thus, the number of dimensionless groups necessary in equation (3) was 
reduced by six/ 

Since the propellant temperature is the lower limit of the ullage mean temperature, 
it is logical to replace all temperatures with temperature differences above the liquid 
propellant temperature (T L ). 

Equation (3) may now be expressed in terms of the following ten dimensionless 
groups (the group numbers have been changed to run from one to ten); 


* l = 


*z = 


T - T 
m B 

T - T t 
o L 


Jg M 2 k (T - T ) 
c w ' o L 

p 3 r 4 


(29) 

(30) 
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7T 3 = 



7T 4 = 


M r 0 T 

M ~~ 
w 


(31) 


(32) 


*5 = 


M k 
w 


P ( C p d ) 

1 ' KITI7 1 TT T 


pw w w 


g M p V. 
c w 1 

“ p 2 r 4 


g M p 
c w 
*7 ~ ,,2 


p “ r 


TTg = 


h r 
a 


TTq = 


T - T 
a L 

T - T t 
o L 


*10 ~ 


M V T 
w L 

V r 2 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


Curve Fit of Dimensionless Equation 

These ten dimensionless groups (equations 29 through 38) can be used to 
correlate the results of tests and computer runs according to equation (3). In most 
cases equation (3) would be written in the following form: 


_ Sy6e4^ T<7 & 

*1 ~ a 7r 2 *3 *4 *5 *7 7 r 8 TTg 7 T 4 o 


(39) 


However, in this case it is necessary to satisfy certain boundary conditions that cannot be 
satisfied by equation (39). The ullage mean temperature at cutoff must remain finite and 
not equal to zero as the ambient heat transfer approaches zero. This boundary condition 
cannot be satisfied by equation (39), but it can be satisfied if the functional dependence on 
7Tg and 7 t 9 is exponential. Also, as the pressur&nt gas distributor Reynolds number 7 r 10 
approaches zero, the heat transfer in the tank approaches free convection. Therefore, 
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the boundary condition of finite, non-zero mean temperature, when tt 10 is zero imposed, 
dictates an exponential functional dependence on ir 10 . Thus, these boundary conditions 
can be satisfied by writing equation (3) in the form 


_ j3 y 6 e 

7Tji — Oij 7T 2 7T 3 ir 4 7T S Tig 


I 4 («“*** *■) (e“ 5,r ‘«) 


(40) 


The coefficients and exponents in equation (40) were evaluated by a curve fit to 
the data from the computer runs and tests. It was found that the data could be correlated 
by this equation if the coefficients a t and in the exponentials were taken as functions 
of 7 t 2 and ir 3 : 


l LO ip 

= “2 n 2 n 3 



(41) 

i (p 

«3 = a 2 H 



(42) 

Equation (40) then becomes 




/3 y 6 e £ 

7Ti = Ol x 7r 2 7T 3 7T 4 TTg 7Tg 

X 

tt 7 

^ e a 2 nf 7 rf 7Tq t r 9 j ^Q! 3 tt? 7r 10 j 

(43) 

where all coefficients and exponents 

are constants. 



It was convenient to divide several groups of equation (43) by powers of ten to obtain 
numbers more easily handled. Therefore, the final equation used in the curve-fit was 



From the curve-fit, the following values were obtained for the coefficients and 


exponents in equation (44): 


= 0. 424 

i = 0. 01416 

a 2 = 0. 00210 

X = 0. 0620 

a 3 = -0. 0292 

co = 0. 415 

13 = -0. 1322 

ip = 1. 174 
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y = -0. 1688 r = 0. 765 

6 = -0. 1146 </> = 0. 1510 

e= 0. 0780 

T m “ T l 

Therefore, since tij — — — — — — , equation (44) becomes 

o L 

™ m -0.1322 -0.1688 -0.1146 

^= 0.424 ^) („) («) 


0.0780 0.01416 0.0620 



This equation is general and is capable of predicting the ullage mean temperature, and 
thus pressurant mass at cutoff within 10 percent for cylindrical tanks with rounded bulk- 
heads. Figure 1 shows total pressurant requirements obtained by various investigators 
for a wide range of tank sizes and system parameters compared with the pressurant 
weights calculated by equation (44). Excellent agreement is obtained over the entire 
range of conditions for hydrogen and oxygen pressurization. * However, the equation is 
limited in its application to conditions of constant ullage pressure, pressurant inlet 
temperature, and ambient heat transfer. The studies indicated that the equation is inac- 
curate at inlet temperatures less than 100° R above the saturation temperature, at ullage 
pressure below propellant saturation, and for very short expulsion times of less than 50 
seconds. The restriction to cylindrical tanks can be removed by proper choice of the 
characteristic tank radius. 


* Evaluating the test results of Reference 1 by this method resulted in a large, but con- 
stant deviation from actual observed pressurant weight. This is probably because the test 
parameters such as heat leak through the vacuum chamber and pressurant inlet tempera- 
ture had to be assumed. Additional information about these tests are required to 
re-evaluate these conditions. 
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Studies have shown that the characteristic tank radius for oblate spheroids, used 
in equation (44) , should be about two-thirds of the maximum tank radius. This assump- 
tion is theoretically justified, because a cylinder having the same volume and surface 
area as an oblate spheroid has a radius equal to 0. 63 times its maximum radius. Further 
test data and analytical studies are necessary to select the characteristic radius for other 
geometries. 

Due to the dimensionless nature of equation (44) , it is not restricted to any 
particular propellant, pressurant, or tank size as indicated in Figure 1. Although an un- 
insulated propellant tank was assumed in the development of this equation, Figure 1 
shows good agreement with test results obtained with vacuum jacketed liquid hydrogen 
tanks. 


To simplify the use of equation (44) , the pressurant and propellant properties 
for the case of liquid oxygen pressurized by oxygen and liquid oxygen pressurized by 
helium were substituted in equation (44) , and the following equations were obtained. 

Since these equations are dimensional they are applicable only to the case indicated. For 
liquid oxygen pressurized by oxygen, 


T - 164 
m 


T - 164 
o 


= 3.33 (T 
' o 


- 164) 


-0.297 


0. 1395 
r 


0.01416 

i 


. -0.0780 0.0762 „ -0.1146 

(C p d ) p 0_ 

pw w w T 




exp | 0.001568 (T q - 164) “°* 304 r °* 895 0,765 (T & - 164) 


■] 


• exp 


. nn „ m , „ . -0. 574 -2.604 f V L 

-120. 9 (T - 164) r 

' "D. 


(45) 


and for liquid oxygen pressurized by helium, 


T - 164 

T “- ~ 64 = 3 - 10 (T o - 164 > 


-0.304 0.1395 „ 0.01416 

r V. 

l 


„ , -0.078 0. 0762 „ -0.1146 

(C p d ) p e 

pw w w T 
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• exp | 0. 0000681 (T q - 164) °' 1650 r ' 0 ’ 895 °- 765 164) 


>] 


exp 


-5.04 (T - 164) -°- 443 r- 2 ’ 604 -± 

° a d 


( 46 ) 


The units of all variables in equations (45) and (46) must be those given in the 
Definition of Symbols. For other combinations of propellant and pressurant, equation 
(44) must be used. After the ullage mean temperature at cutoff is calculated, and using 
equation (44) , (45) , or (46) , if applicable, the total mass of required pressurant gas 
may be calculated from equation (1). 


In designing a launch or space vehicle pressurization system, vehicle parameters 
Such as tank volume, engine flowrate, tank material, etc. , determined by vehicle mission 
profile, are fixed input values. However, there are various controllable parameters in a 
pressurization system that can be used to optimize the system without affecting basic ve- 
hicle characteristics. The relative significance of various parameters on pressurant 
requirements has, therefore, been investigated under another study program. The re- 
sults of these studies exerpted from reference 1 are presented in Figure 2. From a 
central origin, representing a reference condition (Saturn V, S-IC Stage) for all param- 
eters, the increase (+Y) and decrease (-Y) , of the ullage mean temperature at cutoff is 
shown as a function of variation of the parameters on the abscissa. The parameters 
were varied over a range expected for vehicle design. Thus, pressurant inlet tempera- 
ture can increase or decrease by a factor of 2 from the reference condition, pressure 
by a factor of 3, tank radius by a factor of 2, expulsion time by a factor 3, etc. It was 
indicated that the pressurant inlet temperature exerts the greatest influence on the ullage 
mean temperature. Diminishing return of this effect did not exist within the range of in- 
vestigation (530°R to 1200°R). The mean temperature increased as the ullage pressure 
was increased and also as the tank radius was increased. Increasing the tank wall thick- 
ness, heat capacity, or density caused decrease in the mean temperature. The pres- 
surant distributor flow (Aq) that controls the gas-to-wall forced convection heat transfer 
coefficient had a significant effect on the mean temperature when Ajj was reduced, but no 
effect at all when flow area was increased. This indicates that the pressurant inlet ve- 
locity for the reference systems was chosen at an optimum point. Figure 2 also indicates 
that helium pressurant must be introduced into a tank at a temperature 1. 1 times higher 
than oxygen pressurant to obtain the same ullage mean temperature. 


CONCLUSIONS 


An equation derived by dimensional analysis provides a reasonably accurate 
method for prediction of pressurant requirements for cylindrical LOX and hydrogen pro- 
pellant containers. This method is advantageous for preliminary design and optimization 
studies where the use of large computer programs becomes, excessive in cost and time. 
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FIGURE 1. COMPARISON OF CALCULATED AND MEASURED PRESSURANT MASS FOR 
VARIOUS TEST PARAMETERS 
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FIGURE 2. THE EFFECTS OF VARIOUS DESIGN PARAMETERS ON THE MEAN TEMPERA' 
TURE AT CUTOFF 
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